Analysis of the loop length distribution for the negative weight percolation problem in 

dimensions d = 2 through 6 
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We consider the negative weight percolation (NWP) problem on hypercubic lattice graphs with 
fully periodic boundary conditions in all relevant dimensions from d = 2 to the upper critical 
dimension d = 6. The problem exhibits edge weights drawn from disorder distributions that allow 
for weights of either sign. We are interested in in the full ensemble of loops with negative weight, 
i.e. non-trivial (system spanning) loops as well as topologically trivial ("small") loops. The NWP 
phenomenon refers to the disorder driven proliferation of system spanning loops of total negative 
weight. While previous studies where focused on the latter loops, we here put under scrutiny the 
ensemble of small loops. Our aim is to characterize -using this extensive and exhaustive numerical 
study- the loop length distribution of the small loops right at and below the critical point of the 
hypercubic setups by means of two independent critical exponents. These can further be related to 
the results of previous finite-size scaling analyses carried out for the system spanning loops. For the 
numerical simulations we employed a mapping of the NWP model to a combinatorial optimization 
problem that can be solved exactly by using sophisticated matching algorithms. This allowed us to 
study here numerically exact very large systems with high statistics. 
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I. INTRODUCTION 

The statistical properties of lattice-path models on 
graphs, equipped with quenched disorder, have experi- 
enced much attention during the last decades. They have 
proven to be useful in order to characterize, e.g., linear 
polymers in disordered/random media [ll-QjVortex loops 
in high-Tc superconductivity at zero field [6|t2| and the 
(i = 3XYmodel0,IIl| , networks of vortex strings found 
after a symmetry-breaking phase transition in field the- 
ories [l?fl3 | , as well as domain wall excitations in disor- 
dered media such as spin glasses [l^ [1^1 and the solid- 
on-solid model The precise computation of these 
paths can often be formulated in terms of a combinato- 
rial optimization problem and hence might allow for the 
application of exact optimization algorithms [18j] devel- 
oped in computer science. 

For an analysis of the statistical properties of these 
lattice path models, geometric observables and scaling 
concepts similar to those developed in percolation theory 
[19l - [2l| have been used conveniently. In the past decades, 
a large number of percolation problems in various con- 
texts have been investigated through numerical simula- 
tions. Among those are problems, where the fundamental 
entities are string-like, similar to the lattice-path models 
mentioned in the beginning, rather than clusters consist- 
ing of occupied nearest neighbor sites as in the case of 
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usual random bond percolation. 

In a sequence of recent articles we have introduced p2j 
and investigated (see below) the negative-weight perco- 
lation (NWP), a problem with subtle differences as com- 
pared to other string-like percolation problems. In the 
most basic NWP setup, one considers a regular lattice 
graph with periodic boundary conditions (BCs), where 
adjacent sites are joined by undirected edges. Weights 
are assigned to the edges, representing quenched random 
variables drawn from a distribution that allows for edge 
weights of either sign. The properties of the weight dis- 
tribution are further controlled by a tunable disorder pa- 
rameter, signified p. For a given realization of the dis- 
order, one then computes a configuration of loops, i.e. 
closed paths on the lattice graph, such that the sum of 
the edge weights that build up the loops is minimal and 
negative. As an additional optimization constraint we 
impose the condition that the loops are not allowed to 
intersect; consequently there is no definition of clusters 
in the NWP model. Regardless of the spacial dimension 
of the underlying (hypercubic) lattice graph, the observ- 
ables are always line-like, i.e. have an intrinsic dimension 
of c? = 1. Nevertheless, the loops may be fractal with 
fractal dimensions d/ > 1, see Ref. (23j . 

The problem of finding these loops numerically can be 
cast into a minimum-weight path (MWP) problem, out- 
lined in sect. |lT] in more detail. A pivotal observation is 
that, as a function of the disorder parameter p, the NWP 
model features a disorder driven, geometric phase tran- 
sition, p2l - [2^ triggered by a vital change of the typical 
loop size (as discussed below in more detail). In this re- 
gard, depending on the precise lattice setup and on the 
value of p, one can identify two different phases: (i) a 
phase where the loops are "small" , meaning that the lin- 
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FIG. 1: Samples of minimum weight configurations of loops 
for a 2D square lattice with side length L = 64 and fully pe- 
riodic boundary conditions. The snapshots relate to different 
values of the disorder parameter p, where (a) p ~ pc, (b) 
p' < pc, and, (c) p" < p' . In the limit of large system sizes 
and above the critical point pa, loops might span the lattice 
along at least one direction as, e.g., the gray loop in (a). For 
small values of p, loops with a comparatively large length 
appear to be suppressed exponentially. 



ear extensions of the loops are small in comparison to the 
system size, see Figs.[ljb-c) (therein, the linear extension 
of a loop refers to its projection onto the independent lat- 
tice axes), and, (ii) a phase where "large" loops exist that 
span the entire lattice, see Fig.[TJa). Regarding these two 
phases and in the limit of large system sizes, there is a 
particular value of the disorder parameter, signified as 
Pc at which system spanning (or "percolating") loops 
appear for the first time. 

Previously, we have investigated the NWP phe- 
nomenon for 2D lattice graphs i22j using finite-size scal- 
ing (FSS) analyses, where we characterized the under- 
lying transition by means of a set of critical exponents. 
Considering different disorder distributions and lattice 
geometries, the exponents where found to be universal in 
2D and clearly distinct from those describing other perco- 
lation phenomena. In a subsequent study we investigated 
the effect of dilution on the critical properties of the 2D 
NWP phenomenon 



24| . Therefore we performed FSS 



analyses to probe critical points along the critical line in 
the disorder-dilution plane that separates domains that 
exhibit or do not exhibit system spanning loops. One 
conclusion of that study was that bond dilution changes 
the universality class of the NWP problem. Further we 
found that, for bond-diluted lattices prepared at the per- 
colation threshold of 2D random percolation and at full 
disorder, the geometric properties of the system span- 
ning loops compare well to those of ordinary self-avoiding 
walks. We performed further simulations for the NWP 
model on hypercubic lattice graphs in dimensions d — 2 
through 7 [23|, where we found evidence for an upper 
critical dimension d„ = 6 of the NWP phenomenon. This 
result was based on monitoring the critical exponents re- 
lated to the NWP transition (one expects them to stay 
fixed for c? > (i„). We also studied numerically as well as 
analytically a variant of the NWP transition on 3-regular 
random graphs (RRGs) , i.e. graphs where each node has 
exactly 3 neighbors and where there is no regular lattice 
structure. Hence, we obtained direct access to the mean- 
field exponents that govern the model for d > d^- We 
obtained excellent agreement between numerical and an- 



alytic results and could provide further support for the 
claim du = 6- 

All of the studies mentioned above where focused on 
the statistical properties of the largest loop (or more pre- 
cisely, the longest loop) for a given realization of the dis- 
order and the critical properties of the NWP model that 
derive from an analysis of these loops in the vicinity of 
the critical point pc- Up to now, limited attention was 
payed to the ensemble of "small" loops that actually com- 
prise the major part of loop segments in the vicinity of 
Pc- As we found earlier, at this critical point the loops 
are rather isolated and well separated from each other, 
resembhng a dilute gas of loops (cf. Fig. [T]). Further, 
the normalized and ensemble-averaged probability mass 
function (pmf) ni of loop s having length £ right at pc was 
studied in Refs. [l^, |23| . It exhibits an algebraic decay 
similar to the distribution of cluster sizes at the critical 
point in ordinary random percolation (l9l [20|. i.e. 



ni{pc) Oil at p = Pc 



(1) 



The numerical values of the decay exponent t (also 
termed "Fisher exponent" ) found for the NWP model in 
d = 2 through 7 are listed in Tab. HI Note that the Fisher 
exponent is only one out of two independent exponents 
that characterize the whole ensemble of loops. 

In the present article, the second critical exponent that 
characterizes the ensemble of small loops is addressed. In 
this regard, the present article discusses the pmf ni{p) as 
function of the disorder parameter p. Consequently, the 
numerical effort to obtain these distributions, in several 
dimensions d = 2, . . . , 6 was much larger, compared to 
the previous studies where the distribution was obtained 
just for p fa Pc- One of our main results is that for 
values p < Pc the pmf appears to scale similar to the 
distribution of cluster sizes in usual percolation [l^, [l^ , 
i.e. 



ni{p)(x£ ^ exp{-rL(p)^} for p < Pc 



(2) 



TABLE I: Critical properties that characterize the NWP phe- 
nomenon in d = 2 . . . 6. From left to right: Lattice dimension 
d, critical point pc, product i/p ■ d/,p of the critical exponents 
Vp and df^p that describe the divergence of a typical length 
scale and the fractal dimension d/,p (measured at pc), respec- 
tively, as well as the length fluctuation exponent 7p. Further, 
the table lists the Fisher exponent r and the loop-length cut- 
off exponent a. Note that the figures in all but the two right 
columns are taken from Ref ^23]. The last column is meant to 
check the scaling relation 7p = (3 — r)/(T. The additional sub- 
script p indicates that these exponents result from an analysis 
of the percolating loops. 



d L 



pc 



(3-r)/a 



2 5f2 0.340(f) 0.53(3) 0.77(7) 2.59(3) 0.53(3) 0.77(f0) 

3 64 0.f273(3) 0.69(2) -0.09(3) 3.07(f) 0.7f(f) -O.fO(f) 

4 2f 0.0640(2) 0.78(3) -0.66(5) 3.55(2) 0.78(2) -0.7f(f) 

5 f2 0.0385(2) 0.86(4) -f.06(7) 3.86(3) 0.88(2) -0.98(f) 

6 6 0.0265(2) f.00(3) -0.99(3) 4.00(2) 0.97(4) -f.03(2) 
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Therein, the exponential factor accounts for the obser- 
vation that below the critical point Pc the proliferation 
of "long" (still non-spanning) loops is suppressed due 
to some finite "loop size cut-off parameter". The lat- 
ter might be captured by means of a scaling parameter 
Tl(/o) which depends on the subtleties of the dis- 
order. Its inverse fo(p) = ^/Thip) relates to a typical 
length scale to which the perimeter of the loops is lim- 
ited at a given value of p and it should not depend on 
the side length L of the system (at least in the limit of 
large system sizes where a loop of, say, length fits well 
into the simulation box). Therefore, loop configurations 
alSHSl small values of p are consistent with a spanning 
probability Pl{p < Pc) ^ in the hmit of large sys- 
tem sizes (L — cx)). As the critical point is approached 
from below the loop size cut-off parameter vanishes, giv- 
ing rise to the purely algebraic decay of n£ observed at 
Pc as in Eq. ([1]), featuring loops with length £ on virtu- 
ally all length scales. (A qualitatively similar observation 
in the context of high-Tc superconductors is referred to 
as "Onsager vortex-loop unbinding" that signals the su- 
perconductor to normal metal transition Further, 
in string theory, the analo g ob servation is referred to as 
"Hagedorn transition" (l^. Il3l|.') The decrease of the pa- 
rameter Tl(p) can be related to a second, independent 
exponent that, in addition to r, serves to characterize 
the ensemble of small loops. The respective critical ex- 
ponent a is defined via 

Tl(p)(x|p-p,|1/- , (3) 

where a might be referred to as "loop-size cut-off" ex- 
ponent (i.e. the critical exponent related to the loop-size 
cut-off parameter Tl) and where p approaches pc from 
below. Similarly, the corresponding lengthscale £o, to 
which the loops are confined, diverges. This implies that 
loops might get arbitrarily long, limited only by the fi- 
nite size of the underlying lattice. One might expect a 
maximal loop length of imax^ L'^^ , where df denotes the 
fractal scaling dimension of the loops. Thus, at pc and in 
the limit i — > oo , the distribution of the loop perimeter 
exhibits an algebraic decay, solely governed by the fisher 
exponent r. Finally, according to scaling theory ^19, 20], 
the scaling relations 

7p = (3-T)/a (4b) 

should hold, relating <t (as measured for the small loops) 
to Vp, df^p and 7p (all measured from the system span- 
ning loops; indicated by the subscript p). These three 
exponents signify the critical exponents that describe the 
divergence of the correlation length, the scaling dimen- 
sion of the loops, and the fiuctuations of the loop order- 
parameter, respectively. 

The remainder of the present article is organized as 
follows. In section |lTl we introduce the model in more 
detail and we outline the algorithm used compute the 
loop configurations. In section Hill we list the results of 



our numerical simulations and in section ITVl we conclude 
with a summary. Note that an extensive summary of this 
paper is available at the papercore database [25j . 

II. MODEL AND ALGORITHM 

In the present article we consider hypercubic lattice 
graphs G = (V, E) with side length L and fully periodic 
boundary conditions (BCs) for all relevant dimensions 
d — 2 through 6. The considered graphs have A^ = |F| = 
L'^ sites i G V and a number of \E\ = zN/2 undirected 
edges {i,j} & E that join adjacent sites i, j G V. Above, 
z signifies the coordination number of the lattice geom- 
etry, where z = 2d. We further assign a weight ujij to 
each € E. These weights represent quenched ran- 

dom variables that introduce disorder to the lattice. Here 
we consider independent identically distributed weights 
which either have just weight one (probability 1-p) or 
are drawn (probability p) from a Gaussian distribution 
with zeor mean and veriance one. Hence, the disorder 
distribution is given by 

P(cj) = pexp(-cjV2)/V2^+(l-p)(5(a;-l), (5) 

that explicitly allows for loops C with a negative total 
weight ^c — ^s^ij^^c^ij- '^'^ support intuition: For any 
nonzero value of the disorder parameter p, a sufficiently 
large lattice will exhibit at least "small" loops that ex- 
hibit a negative weight, see Fig. [UJc). If the disorder 
parameter is large enough, system spanning loops with 
negative weight will exist, see Figs.[lja). 

The NWP problem then reads as follows: Given a re- 
alization of the disorder for a hypercubic lattice graph 
G, determine a set C of loops such that the configura- 
tion energy, defined as the sum of all the loop-weights 
E = "YIiCgc '^■C' minimized. As further optimization con- 
straint, the loops are not allowed to intersect. Note that 
due to the "energy minimization principle" of the opti- 
mization procedure, the weight of an individual loop is 
necessarily smaller than zero. The configuration energy 
£ is the quantity subject to optimization and the result 
of the optimization procedure is a set of loops C, obtained 
using an appropriate transformation of the original graph 
[26| . For the transformed graphs, minimum-weight per- 
fect matchings (MWPMs) [27 23| are calculated, that 
serve to identify the loops for a given realization of the 
disorder. Since exact MWPMs can be obtained in poly- 
nomial time, this procedure allows for an efficient imple- 
mentation [3(1 of the simulation algorithms. Here, we 
give a brief description of the algorithmic procedure that 
yields a minimum-weight set of loops for a given real- 
ization of the disorder. Fig. [2] illustrates the three basic 
steps, detailed below: 

(1) each edge, joining adjacent sites on the original 
graph G, is replaced by a path of 3 edges. Therefore, 
2 "additional" sites have to be introduced for each edge 
in E. Therein, one of the two edges connecting an ad- 
ditional site to an original site gets the same weight as 
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FIG. 2: Illustration of the algorithmic procedure: (a) original 
lattice G with edge weights. For clarity, a bimodal distri- 
bution that yields edge-weights ±1 is considered. Further, 
only negative edge-weights are shown. Unlabeled edges have 
weight +1. (b) auxiliary graph Ga with proper weight as- 
signment. Black edges carry the same weight as the respec- 
tive edge in the original graph and gray edges carry zero 
weight, (c) minimum- weight perfect matching (MWPM) M: 
bold edges are matched and dashed edges are unmatched, 
and (d) loop configuration (bold edges) that corresponds to 
the MWPM depicted in (c). 



the corresponding edge in G. The remaining two edges 
get zero weight. The original sites i £ V are then "dupli- 
cated", i.e. i — > ii, 12, along with all their incident edges 
and the corresponding weights. For each of these pairs 
of duplicated sites, one additional edge {ii, 12} with zero 
weight is added that connects the two sites ii and 12- 
The resulting auxiliary graph Ga = (Va, Ea) is shown in 
Fig. EJb), where additional sites appear as squares and 
duplicated sites as circles. Fig. [2][b) also illustrates the 
weight assignment on the transformed graph Ga- Note 
that while the original graph (Fig. Hl^a)) is symmetric, 
the transformed graph (Fig. Eljb)) is not. This is due to 
the details of the mapping procedure and the particu- 
lar weight assignment we have chosen. A more extensive 
description of the mapping can be found in . 

(2) a MWPM on the auxiliary graph is determined 
via exact combinatorial optimization algorithms [31]. A 
MWPM is a minimum- weighted subset M of Ea, such 
that each site contained in Va is met by precisely one 
edge in M. This is illustrated in Fig. EJc), where the 
solid edges represent M for the given weight assignment. 
The dashed edges are not matched. Due to construction, 
the auxiliary graph consists of an even number of sites 
and the transformation procedure described in step (1) 
guarantees that a perfect matching exists. Note that a 



MWPM can be computed in polynomial time as a func- 
tion of the number of sites, hence large systems with 
hundreds of thousands of sites are feasible. 

(3) finally it is possible to find a relation between the 
matched edges M on Ga and a configuration of negative- 
weighted loops C on G by tracing back the steps of the 
transformation (1). As regards this, note that each edge 
contained in M that connects an additional site (square) 
to a duphcated site (circle) corresponds to an edge on 
G that is part of a loop, see Fig. ^d). Note that, by 
construction of the auxiliary graph, for each site ii or 
i2 matched in this way, the corresponding "twin" site 
i2/ii must be matched to an additional site as well. This 
guarantees that wherever a path "enters" a site of the 
original graph, the paths also "leaves" the site, corre- 
sponding to the defining condition of loops. All the edges 
in AI that connect like sites (i.e. duplicated-duplicated, 
or additional-additional) carry zero weight and do not 
contribute to a loop on G. Once all loop segments are 
found, a depth-first search [H, can be used to identify 
the loop set C and to determine the geometric proper- 
ties of the individual loops. Here, the exemplary weight 
assignment illustrated in Fig. [SJa) yields 2 loops, i.e. 
C — {£1, £2}, with weights tuci = '-0C2 = ~4 and lengths 
£1 — J2{i j}eCi 1 = ^2 = 4. Hence, the configurational 
energy reads £ = — 8. 

The result of the calculation is a collection C of loops 
such that the total loop weight, and consequently the 
configuration energy £, is minimized. Hence, one obtains 
a global collective optimum of the system. Obviously, all 
loops that contribute to C possess a negative weight. Also 
note that the choice of the weight assignment in step (1) 
is not unique, i.e. there are different possibilities to choose 
a weight assignment that all result in equivalent sets of 
matched edges on the transformed lattice, corresponding 
to the minimum-weight collection of loops on the origi- 
nal lattice. Some of these weight assignments result in a 
more symmetric transformed graph, see e.g. (26j . How- 
ever, this is only a technical issue that does not affect 
the resulting loop configuration. Finally, for the purpose 
of illustration, a small 2D lattice graph with free BCs 
was chosen intentionally. The algorithmic procedure ex- 
tends to higher dimensions and fully periodic BCs in a 
straight-forward manner. 

In the remainder of the article, we will use the pro- 
cedure outlined above in order to study the pmf of loop 
lengths in the NWP model on hypercubic lattice graphs 
in d = 2 through 6. 



III. RESULTS 

Within our extensive numerical studies, we performed 
exact NWP loop calculations for dimensions d — 2 
through 6 for various values p < Pc while averaging over 
many realizations of the disorder. Details are given in 
Tab. nil 

In order to get a grip on the loop-size cut-off parameter 
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a for a particular hypercubic lattice setup of dimension d, 
the pmf n^{p) of the loop length needs to be obtained for 
different values of the disorder parameter p< Pc- Then, 
a best fit to the form ng{p) — no£~'^ exp{— rL(p)^} might 
be used to obtain the three fit parameters no, r, and 
Ti,{p) for different values of p. Finally, the sequence of 
fit parameters T'l(p) might be analyzed to yield the ex- 
ponent <J according to Eq. ([3]). 

However, a different procedure appears to be more ap- 
pealing: from previous simulations in dimensions d = 2 
through 7, reported in Refs. 2^, we found that the 
pmf ni{pc) exhibits an algebraic decay governed by the 
exponent r. For the largest lattice graphs simulated for 
the various dimensions d, we obtain the numerical esti- 
mates listed in Tab. U] For the corresponding data analy- 
ses, very small loops have to be neglected since they are 
affected by the granularity of the lattice and very large 
loops have to be withdrawn since they are affected by 
the lattice boundaries. Once the exponent t for a given 
dimension d is obtained, it can be utilized in the analysis 
of the loop length pmf at p < pc to limit the number of 
fit parameters to only two (i.e. tiq and Tl), allowing for 
a more precise estimate of the individual values of Tl (p) . 

So as to get a grip on the loop size cut-off parameter cr, 
the loop perimeter distribution was obtained for different 
values of the disorder parameter p < Pc and the data 
was fitted using a function as given in Eq. ([5]), with an 
additional normalization factor. E.g., in d = 2 we use 
T — 2.59 fixed, see Tab. |T] and Fig. [31 wherein the fit 
intervals for individual values of p where restricted to 
the range [Anin, ^max]- We further fixed ^min = 10 and 
^max = 60 at p = 0.24 (the upper bound shifting up to 
^max = 150 at Pc). The resulting values for the cut-off 
parameter were then analyzed using a three parameter 
fit to functions of the form of Eq. i.e. Ti^{p) = A\p — 
^cff|i/cr^ Therein the amplitudes A are not of interest 
and the effective critical points p^J^ can be expected to 
differ slightly from the asymptotic critical points listed in 
Tab. U (see discussion below). The resulting loop- length 
cut-off parameters a are also listed in Tab. HI 

For d = 2 systems of side length L = 512 the anal- 
ysis yields a = 0.53(3) and pf = 0.344(2), see Fig. H 
(at L = 256 we find pf = 0.346(2) and a = 0.53(2)). 
For all values of p considered, the data curves of Tl at 



TABLE 11: Simulation parameters: we performed our study 
for Tip values of the disorder parameter p in intervals [pi , p-z\ 
and for a number ur of realizations, for the different dimen- 
sions d and system sizes L. 
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FIG. 3: Results for d — 2 square lattice graphs. The main 
plot shows the probability mass functions ni{p) of the loop 
perimeter £ for different values of the disorder parameter p 
for square systems of side length L = 512. The data curves 
illustrate the suppression of loops with a perimeter larger than 
some cut-off length scale £o, related to a finite "loop size cut- 
off parameter" Ti,(p) for p < pc (see text). The data curves 
fit well to functions of the form of Eq. 0. Right at pc the 
data curve exhibits a pure algebraic decay according to Eq. 
([T]). The inset compares the resulting values of the cut-off 
parameter Ti,{p) for systems of side length L — 512 and L = 
256, where the solid line indicates a best fit to the L — 512 
data using a function with four free parameters as explained 
in the text. 

L = 256 and L = 512 compare well as shown in the 
inset of Fig. [3l This finding is further consistent with 
the usual scaling relation l/{vpdf^p) = a that relates a 
to i/p and df,p, where the latter two exponents signify 
the critical exponent that describes the divergence of 
the correlation length and the scaling dimension of the 
loops at the critical point, respectively. From the re- 
spective values previously obtained '25] one readily finds 
l/(i'pC?/_p) — 0.53(3). Only the location of the (effec- 
tive) critical point pf — 0.344(2), as estimated from 
the complete ensemble of loops at the particular sys- 
tem size L = 512, differs slightly from the one obtained 
from the sca ling analysis of the percolating loops, i.e. 
Pc = 0.340(1) [221 (a similar effect was found in the anal- 
ysis of a d = 3 vortex loop network in Ref. [l3|)- Bear- 
ing in mind that the latter value represents an extrap- 
olation to the thermodynamic limit, the aforementioned 
difference is likely due to finite-size effects. In this re- 
gard it does not come as a surprise that an analysis 
of the loop size cut-off parameter according to Eq. Q 
for fixed pf = 0.340 (instead of pf = 0.344) yields 
the exponent a = 0.61(1) which significantly overesti- 
mates the result obtained from the three parameter fit 
reported above. Finally, a four parameter fit according 
toTi^ip) = Tl + A\p-pf\^^'' at L = 512 (see inset of Fig. 
^ results in the estimates = 0.001(1), pf = 0.348(6) 
and (7 = 0.55(3), with the latter two fit parameters in 
agreement with the ones found above and iii agree- 
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FIG. 4: Scaling analysis of the loop size cut-off parameter 
Ti^{p) in dimensions d — 2 through 6. The solid lines indicate 
fits to functions of the form Tl{p) = A\p ~ pfS^-l". Therein 
the amplitudes A are not of interest and the effective critical 
points pj;^ can be expected to differ slightly, but within sta- 
tistical error bars, from the asymptotic critical points listed 
in Tab. [H The resulting loop-length cut-off parameters a are 
also listed in Tab. U 

ment with zero as one would naively expect. 

So as to facilitate a qualitative comparison, considering 
d = 3 hypercubic lattices with L = 64 we obtained the ef- 
fective critical point p^^ — 0.1278(1) which slightly over- 
estimates the asymptotic value of pc — 0.1273(3) (similar 
to what we observed above in 2Z3), see Fig. The scal- 
ing exponent a — 0.71(1) compares well to the product 

■ ^/,p ~ 0.69(2) for the respective dimension. In the 
analysis of the data for d > 3 we obtained most satisfac- 
tory fits by fixing the parameters to their expected 
asymptotic values pc listed in Tab. HI The analysis of the 
loop size cut-off parameter in dimensions d = 2 through 
6 are shown in Fig. 3] and summarized in Tab. ID 

In higher dimensions we also checked that the loop size 
cut-off parameter T'l(p) at a given value of p is practi- 
cally independent of the system size (however, for larger 
system sizes the loop yield is bigger and hence the statis- 
tics get more reliable). E.g., for d = A hypercubic lattices 
at p = 0.026, i.e. with some distance to the critical point 
Pc = 0.0640(2), we obtained Tl(0.026) = 0.43(5) at L 
16 and Tl(0.026) = 0.45(4) at L = 21. Close to the criti- 
cal point at p = 0.058 we further find rL(0.058) = 0.35(1) 
at L = 16 and rL(0.058) = 0.036(1) at L = 21. Further, 
we observe that in any dimension considered, the scaling 
relation 7p = (3 — T)/a appears to be satisfied within 
errorbars. 



IV. CONCLUSIONS 

In the presented analysis of the NWP model, we 
performed numerical simulations on hypercubic lattice 
graphs for all dimensions relevant for the model, i.e. d=2 



through 6. The aim of the study was to characterize the 
ensemble of small loops in the NWP model using two 
independent critical exponents: the Fisher exponent t 
(which was already known from previous studies, see Ref. 
,23j), and the loop- length cut-off exponent a. Both expo- 
nents can be determined by means of an analysis of the 
probability mass function ni(p) measuring the distribu- 
tion of loop lengths i, considering a sequence of different 
values p close to but below the critical point. This implies 
a huge numerical effort, since we had to study in different 
dimensions large systems, for several values of the disor- 
der parameter, while averaging over many realizations 
of the disorder. For the numerical simulations we used 
a mapping of the NWP model to a combinatorial opti- 
mization problem that allows to obtain configurations of 
minimum weight loops via exact algorithms. Note that 
due to the small side lengths of the lattice graphs that 
are accessible in high dimensions, the data analysis is no- 
toriously difficult at large values of d. However, we find 
the results regarding the exponent a consistent with the 
scaling relations Eqs. (j4a|) and (|4bp for any dimension 
considered (see Tab.[l|. Thus, via this extensive numeri- 
cal study we have completed a comprehensive description 
of the static behavior of the NWP model in all relevant 
dimensions d = 2, . . . , 6. 

In particular, many of the d — S loop models stud- 
ied in the literature report on values for the critical 
exponents r and a that are close to r = 3.07(1) and 
cr = 0.71(1) (1/cr = 1.41(2)) found here. E.g., Ref. 
obtains t — 2.4(1) (however, in a previous study they 
report t = 3, see Ref. @) and l/cr = 1.45(5) (in that 
study the latter quantity was called 7) for the d = 3 
uniformly frustrated XY model as well as for the lattice 
Ginzburg-Landau model in a frozen gauge approxima- 
tion, and Ref. |81] yields r = 2.8(1) and a = 0.6(1) for the 
strongly screened vortex glass model. Note that in d = 3 
all the critical exponents of the NWP problem appear to 
be quite close to those that describe the strongly screened 
vortex glass model analyzed in Ref. Q. As regards this, 
it appears to be tempting to conclude that in d = 3 both 
models are in the same universality class. 
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